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Abstract. We study valuation of swing options on commodity markets when the commodity 
prices are driven by multiple factors. The factors are modeled as diffusion processes driven by a 
multidimensional Levy process. We set up a valuation model in terms of a dynamic programming 
problem where the option can be exercised continuously in time. Here, the number of swing 
rights is given by a total volume constraint. We analyze some general properties of the model 
and study the solution by analyzing the associated HJB-equation. Furthermore, we discuss 
the issues caused by the multi-dimensionality of the commodity price model. The results are 
illustrated numerically with three explicit examples. 

1. Introduction 

The purpose of this paper is to propose and analyze a model for valuation of a swing option, 
see, e.g. [B], written on multiple commodities when the commodity spot prices are driven by 
multiple, potentially non-Gaussian factors. More precisely, the model is formulated as a dynamic 
programming problem in continuous time. The holder of the option is contracted an amount of 
a given commodity that can be purchased for a fixed price during the lifetime of the contract. 
The purchases can be done (that is, the option can be exercised) continuously in time such that 
contracted rate constraints are fulfilled. This form of contract originates from electricity markets, 
where they are called flexible load contracts, see, e.g. [51 [TH]. However, this model setting can also 
fit a traditional swing option with a high number of swing rights and possible exercise times. For 
example, we can think of a situation in an electricity market where contract is written for a year 
and holder can exercise on the hour- ahead market. This results into over 8000 possible exercise 
times, which makes, in particular, Monte-Carlo methods virtually intractable. 

During the recent years, there has been a lot of activity on analysis of swing options. Being 
essentially a multi-strike American or Bermudan option, a natural way to approach swing options is 
via an optimal multiple stopping problem. In the recent papers [3 [2], the theory of optimal multiple 
stopping is developed in continuous time using sophisticated martingale theory. To compute option 
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prices numerically, they develop appropriate Monte-Carlo methodology. Other methodology for 
swing option pricing includes forests of trees [T31 [131 HOI ttH] or stochastic meshes P2] , multi-stage 
stochastic optimization 12 , (quasi-) variational inequalities [HI US] and PDE approaches [Tni^lTS]. 
Fundamentally, all of these methods are based on the dynamic programming principle. 

As the main contribution of this paper, we develop a valuation model for multi-commodity 
swing options inspired by [4J. In [^, the valuation problem was studied in the case of a single 
commodity driven by a one-factor Gaussian price process. In this paper, we generalize the results 
of [4] to cover multiple contracted commodities with prices driven by multiple factors. From 
applications point of view, this is an important generalization, since there is a substantial body of 
literature supporting the usage of multi-factor models for commodity prices. Moreover, we allow 
also for non-Gaussian factors, which are favored, for example, in electricity price models, see, e.g. 
[21 [13] . We model the factors as a multi-dimensional Levy diffusion and the underlying commodity 
prices are obtained by a linear mapping of the factors. This makes our model more tractable yet 
keeping it still very flexible as it allows us to take, for example, heat rates and spreads into account 
in a natural way. Our study is also related to [15;, where a similar model is used to study hedging 
of swing options. We also refer to 17J, where swing option pricing is considered under a non- 
Gaussian multi-factor price model. However, the analysis of [I7j is restricted to a modification 
of the so-called Deng model (see [S]), which is a particular mean-reverting model. In our paper, 
we set up and analyze a class of models where the underlying factor prices follow a general Levy 
diffusion. The existing mathematical literature on swing options is mostly concerned with the 
pricing of a swing option. In addition to pricing, we also address the question of how to exercise 
a swing option optimally. From the analytical point of view, we identify using the HJB-equation 
an optimal exercise policy and characterize it in an intuitive way using the notion of marginal lost 
option value. We also present a numerical analysis of the problem including a numerical scheme 
based on the finite difference method. 

The reminder of the paper is organized as follows. In Section 2 we propose our model for the 
valuation of swing options. In Section 3 we analyze some general properties of the value function. 
Section 4 is devoted to the derivation of necessary and sufficient conditions for a function to 
coincide with the value function. We illustrate our results with explicit examples in Section 5, 
which are solved numerically in Section 6. Finally, we conclude in Section 7. 
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2. The valuation model 

2.1. The price dynamics. As we mentioned in the introduction, the prices of the commodities 
are driven by multiple factors. Throughout the study, the number of commodities is m and the 
number of driving factors is n. The factor dynamics X are modeled by an n-dimensional Levy 
diffusion. To make a precise statement, let (17, J^, F,P) be a complete filtered probability space 
satisfying the usual conditions, where F = {J-}t>o is the filtration generated by X. We assume 
that the factor process X are given as a strongly unique solution of the Ito equation 

(2.1) dX{t) = a{t,X{t))dt + B{t,X{t))dWit)+ [ ^{t,X{t),ON{dt,dO, 

where W = {W^, . . . , W'^'') is an nf,-dimensional, potentially correlated, Brownian motion satisfy- 
ing d{Wi, Wl) = pi-jdt with pij e [-1, 1] for ah Furthermore, N = (iV\ . . . , A^"') denotes an 
rt(-dimensional Poisson random measure with Levy measure v given by the independent Poisson 
processes 7y^, . . . Here, j^({0}) is the unit measure concentrated on zero and it is finite. The 
coefficients a : [0,T] x M" W\ a : [0,r] x M" ^ K"xn, = . [q^j-] x M" x M"' — > M"><"' 
are assumed to be sufficiently well behaving Lipschitz continuous functions to guarantee that the 



Ito equation (2.1) has a unique strong solution - see [T], p. 365 - 366. The motivation to model 
the randomness using Brownian and finite activity jump noise comes from electricity prices. In 
this framework, the jump process models the spiky behavior in the prices whereas the Brownian 
motion takes care of the small fiuctuations. 

Using the factor dynamics X, we define the m-dimensional price process t ^ Pt :— P{Xt) via 
the linear transformation 

(2.2) P{x) Bx, 

where x G M" and i? is a constant m x n matrix with rank[B) — m < n. In other words, there 
exists constants {bij) such that P*(a;) — X]j=i ^ij^^ foi' all i = 1, . . . , to, that is, the commodity 
prices are linear combinations of the driving factors. The component P' models the time evolution 
of the price of the ith commodity and this price is driven by the n factors, i.e. the n-dimensional 



Levy diffusion X given as the solution of the Ito equation (2.1 1. Since the price is linear as a 



function of the factors it is easy to change the model into a price model for spreads. Furthermore, 



the matrix B in (2.2) can be interpreted as a constant weight between the different factors X 



affecting the price. That allow us to take, for example, heat rates into account in our model. 
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In the definition of the factor dynamics, we assumed that the jump-diffusion X and the driving 
Brownian motion and Levy process have all different dimensions. For notational convenience, we 
assume in what follows that these dimensions are the same, i.e. n — ni, — ni. We point out that 
the following analysis holds with obvious modifications also in the case where these dimensions 
are different. 



2.2. The valuation model. The swing option written on the price process P = P{X) gives 
the right to purchase the given amount M of the commodities i over the time period [0,r]. In 
addition to the global constraint M, the purchases are also subject to a local constraint u which 
corresponds to the maximal number of swing rights that can be exercised on a given time. Since 
the swing option can be exercised in continuous time, the local constraint is the maximum rate at 
which the option can be exercised. To formalize this, let W = Um^u^ be the set of F-measurable, 
real- valued processes ~ satisfying the constraints 

u\ e [0,?!*], / ulds < M\ 
Jo 

for alH = 1, . . . , TO and t G [0, T]. Here, the elements u = (m*) £ 
M™-valued process Z defined as 

(2.3) 



and M = (A/') € 



The 



Zl = I Usds, 



where i = l,...,m, keeps track of the amount purchased of commodity i up to time t. In 
what follows, we call Z the total volume and denote the product as U. The integral 



representation for in (2.31 is well defined due to the local constraint. 

Denote the set S [0, T] x 0" ^[0, Mi] x M" and define the affinc function A 



as 



A{x) ^Qx + K, 



where Q is an to x to matrix and K G M'" . Define the expected present value of the total exercise 
payoff J : S X U ^ R given by the rate u £ hi from time t up to the terminal time T (or, the 
performance functional of u) as 



(2.4) 



E 



1=1 



Zt — Xt — X 



where r > is the constant discount factor. We point out that function J is defined explicitly 
as a function of the factors X. This corresponds to that the holder of the contract observes the 
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underlying factor and bases her exercise decisions of this information. Furthermore, we remark 
that this framework covers essentiahy call- and put-like payoffs, where the strike prices are given 
by the constant vector K. Now, the value function F : 5 — >■ M is defined as 

(2.5) V{t, z,p) — swp J{t, z,p,u). 

ueu 

We denote an optimal rate as u* . 



We make some remarks on the valuation problem (2.5 1. The dimension of the decision variable 
u is the same as the dimension of the price. That is, we can exercise the option for each price 
component, which corresponds to different commodities, with a different decision variable. Fur- 
thermore, we defined the function A such that it takes values in M™. This is done for notational 
convenience. Suppose that we have an m-dimensional price process but the decision variable u is 
/c-dimensional with k < m. This corresponds to the case where m commodities are bundled into 
k baskets and the holder can exercise the option on the baskets. Formally this is done by defining 
the affine function as A : M™ R'^. This will not affect the form of the value function. 



3. Some General Properties 



In this section we study some general properties of the valuation problem (2.5 1. We split the 
analysis in two cases, depending on whether Af ' > u^T or AP < u^T for a given commodity 
i. In the latter case, the limit AP imposes an effective constraint on the usage of the option in 
the sense that the amount A'P is dominated by the amount that can be purchased if the option 
is exercised on full rate over the entire time horizon. This case, i.e. the case when an effective 
volume constraint is present, is the interesting one from the practical point of view. It is also 
substantially more difficult to analyze mathematically as we will see later. Before considering this 
case, we study the complementary case when the effective volume constraint is absent. This will 
give us a point of reference in the other case. 



3.1. Without an effective volume constraint. We consider first the case where KP > u'T for 
a given commodity i. The total volume constraint for the commodity i is now superfluous, since it 
is possible for the holder to exercise the option at full rate throughout the lifetime of the contract. 
In the absence of an effective volume constraint for the commodity i, an optimal exercise rule is 
given by the next proposition. 
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Proposition 3.1. Assume that AP > u^T for a given commodity i. Then an optimal exercise 
rate u** for the commodity i reads as 

u\ tfA'{Pt)>0, 

for allte [0,T]. 

Proof. Let u eU and t G [0,T]. First, we observe that u** GW. Furthermore, we find that 



E 



Zt — z, Xt — X 



E 



(3.1) 



E 



E 



Zt^z, Xt ^ X 



Zt — z, Xt — X 



1=1 

e-r(^-t) A\Pt)t^,l^AHP,)<o}ds 
e~ris-t)A^Pt)umA.^P^)^o}ds 

in 

e''-i''')^A\Pt)u\ds 
1=1 

e~''^'~''>A\Pt)ul*ds 

/ e-'-^'-'^Y^A^iPtHds 
1=1 

Now, take supremum over ah u' on the left hand side and supremum over u\ I i, on the right 



< E 



E 



Zt = z, Xt = X 



Zt — z, Xt — X 



hand side of (3.1). Since the functional J is linear in m, the same inequality still holds and. 



consequently, the conclusion follows. 



□ 



Proposition |3 . 1| states that in the absence of an effective volume constraint for commodity i, it 
is optimal to exercise the option whenever the payoff A*(Pt) is positive, i.e. when {QPtY > —K^. 
This is a natural result, since the holder does not have to worry of running out of the option over 
the planning horizon. Furthermore, since u is a constant vector, we find using Proposition |3.1 
that the value function does not depend on z* in the absence of an effective volume constraint for 
commodity i. This yields the following corollary. 



Corollary 3.2. In the absence of an effective volume constraint for a given commodity i, the 
marginal value ^{t,z,p) = 0. 
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Corollary [3^2] is also a very natural result. Indeed, if the holder uses the option on a commodity 
with no effective volume constraint, the option will not lose value. 

To close the subsection, we discuss how the dimension of the range of the function A affects the 



value given by (2.5). For simplicity, assume that there is no effective volume constraint for any of 



the commodities and that the function A : M™ is of the form 

(3.2) A{x) = diag(gi , . . . , q^n) ■ x + K, 

for K E M™ . Using Proposition |3.1[ we know that the optimal exercise rule for the valuation 



problem specified by the payoff structure (3.2) is 



(3.3) uf 



ifA'(Pt)>0 
ifA'(Pt)<0, 



for all t £ [0,T] and I = 1, . . . ,m. Formally, we can decrease the dimension of the range of A 
from TO, for example, as follows. Take m' < m and define the to' x to- matrix Q — (qij) such that 
each Qi occurs only once and on exactly one column of Q and the other elements are zero. In 
financial terms, this means that the commodities i are bundled into to' pairwise disjoint baskets 
with weights q,;. Then the option gives exercise rights on each of these baskets with separate 
exercise rates. Now, let the function A : M™ M™' be A{x) = Q-x + k with K e M™' such that 

rn 

(3.4) k^^Y.^!^ 

i=i 



for all i — 1, . . . ,to'. Using the same reasoning as in Proposition 3.1 we find that the optimal 



exercise rule for the valuation problem (2.5) given by A is 



(3.5) u 



v} ifi'(Pt)>0 
ifi'(Pt)<0, 



for all t e [0, T] and I = 1, . . . , m'. Denote the value for TO-diniensional (TO'-dimensional) problem 
as V {V). Furthermore, denote the TO'-dimensional total volume variable as Z and assume that 
all maximal exercise rates coincide: u' = = u for all I = 1, . . . , to and I' = 1, . . . , to'. Then, 
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due to the structure of matrix Q, we find using (3.4 1 that 



(3.6) 



Vit,z,p)=-E 



> E 



where the events 



-r(s-t) 



1=1 



Zt — z, Xt 



Vit,z,p), 



in' m 

j=i I j=i 



\ 



ds 



J 



Zt ^ z, Xt = X 



and the cumulative variable Z is defined analogously to (3.4). Summarizing, we have shown that 
by bundling commodities i into mutually disjoint baskets and, thus, reducing the dimension of the 
exercise rate process u, we lower the value of the option. This is, again, a natural result, since 
the bundling of commodities lowers flexibility of option contract in the sense that the holder must 
exercise the option at the same rate for all commodities in the same basket. This is in contrast to 
the case with separate commodities, where the exercise rates can be decided individually for each 
commodity. 



3.2. With an effective volume constraint. In this section, we consider the case where Mi < 
UiT, in other words, the case when the total volume constraint is less than the maximal amount 
of commodity that can be acquired over the lifetime of the option. From the practical point of 
view, this is the interesting case. It is also substantially more difficult to analyze, since in this 



case we cannot find an optimal exercise policy explicitly as in Proposition 3.1 Instead we find the 
value function as the solution to the HJB-cquation and an optimal exercise policy is obtained as 
a biproduct. 



Our first task is to write the conditional expectation in (2.5) such that it depends explicitly 
on Z. This will be helpful in the later analysis. To this end, define the process Y a.s Yt = 
e"''* Elli A\P{^t))Zi. Then the Ito formula yields 



(3.7) dYt = -re-^' E A\P{Xt))Zldt 



1=1 

m m 



i=i 1=1 j=i ^ 
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+ e-''*^yl'(P(Xt))u^di 



1=1 



1 = 1 



d{W\W'') 



Note that since A is affine and P is linear, we have 



(3.8) 



A\P{X,- + 7^=) (t, X,- , e))) - A'(P(X,- )) 

m 
m 

= E9'-^'(^^'H*'^*-'^'))' 



(3.9) 



(3.10) 



9^A'(P(Xt)) 
dx'^dx^ 



= 0, 



5A'(P(X0) aP^(X,) . V- , 



where gy, 6ji and cu are constants for all i,j,L Substitution of (3.81, (3.91 and (3.10) into (3.7| 
yields 



n m 



dYt = -re-^* A\P{Xt))Z\dt + ^ e"''* ^ Z^Q,(a,(t, Xt)dt + a,(i, Xt)dVFt) 

m 

+ e-''*^A'(P(Xt))Mtdt 

n mm 

(3.11) +^ / e-'*5]z,'^ffi„P''(7('=)(i,X,-,C'=))iV'=(di,dC'). 



fc=i • 



;=i v=i 



Since 



(3.12) 



N(dt, di) ^ N{dt, rfO + Kd^dt, 



where v is the Levy measure, we find that (3.11 1 can be written as 



n m 



dYt = -re-^* A'(P{Xt))Z\dt + ^ e"''* ^ Z^q,(a,(t, X,)dt + a,(i, Xt)dT^t) 



1=1 



i=l 1=1 



10 
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1=1 

n m m 

By integrating this from t to T, we obtain 

m m 

m n m 

n » mm 

+ / E / e-'-^E^^E*'^'(7*'H^.^.-.e'))^'(rf^.rfe')- 



1=1 

T 



ds 



(3.13) 



k=l ' 



1=1 v = l 



Consider first the Brownian integral e~^^ ^^^-^^^]^-^ Zlciiai{s, Xs)dWs. Each of the integrands 
is of the form e^^'^ Zlcnaij^s, Xs)dW^ . By definition of Zl, we know that < Z^ < uH. Since Zl 
is nondecreasing, it foUows that (Zlj.)'^ < {uHY < (w'P)^. Hence, 



E 



e-^^^{Z'^^clal{s,Xs)ds 



Zn = Z, Xc)= X 



Zn = Z, Xn ^ X 



< oo. 



/ '7rjis,Xs)ds 

Jo 

Using a martingale representation theorem, see, e.g. [T], Thrm. 5.3.6, we conclude that 

/ e~^'Y.Y.^'sCu(J^{s,Xs)dW, 
is a martingale with respect to F. Using the same argument, we find that the process 

l-t ri 



/ E / e-'-^^Z^^g;.P^'(7('=Hs,X,-,e'=))iV'=(rfs,de'=) 



is a also a martingale with respect to F. Consequently, the conditional expectation with respect 



to Ft is zero for the last two terms in (3.13). 



SWING OPTIONS IN COMMODITY MARKETS 



11 



By multiplying (3.131 with e''* on both sides, substituting into (2.51 and using the martingale 
properties, we find 



V(t,z,p) = supE 



g-r(T-t) J2 A\P{XT))Zir - J2 A\P{Xt))Z\ 



1=1 



1=1 



1=1 



i=l 1=1 



fe=l'^'^ 1=1 v=l J 



(3.14) 



= 1 v=l 



Zt ^ z, Xt = X 



Using that the measurability of J2i=i ^ {P{^t))Zi, we can express the value function (2.5) as 



F(t,z,p) = -E^'b)^' + supE 



1=1 

T / 



ueu 



-r(T-t) 



J2A\PiXT))Z 



1 = 1 



1=1 

m m 



i=l 1 = 1 



IL n lib lib \ 

fc=l'^^ 1=1 v = l ) 



Zt = z, Xt ^ X 



1=1 v = l 

We now have an explicit dependence on Z in the value function, which will be useful in the proof 
of the following proposition. We point out that we can assume that we have an effective volume 
constraint in all commodities i, since the complementary case is already covered by Proposition 



EH 



Proposition 3.3. In the presence of an effective volume constraint, i.e. when M* < u*T, the 
marginal value Vzi{t,z,p) < for all i. 



Proof. Let = (wf , . . . G be processes giving rise to admissible exercise policies Z^ = 
{Z} +£,..., +e) at time t. Let u £U\)e the processes giving rise to admissible exercise policies 
Zt at time t. Since the exercise policies Zf arising from are admissible and must satisfy the 
effective volume constraint we have that lA^ C lA. Also, for an arbitrary admissible s i— > Z| on 
(t, r], define an associated Z as 

(3.15) Zs - Zl - e, 

With this in mind, we proceed by expressing the marginal value as 
-A^{p) 



for all s e (t,T\. 
dV{t,z,p) 



12 
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lim - 

£^■0 S 





( 


< supE 


e-r(T-t) 




V 



Ai{P{XT))Z^.{T) + Y,A\P{Xt))Z\T) 



1=1 



i: 



1=1 



E 

fe=i 



i=l 



-r(s-t) 



1=1 
l¥^3 



( 



\ 



Z]{s)^Y.^\s) 



\ 



1=1 



I 



^Y.^uP''{i^^\s,x,-,e))^\de) 

v=l 

- supE 

u 

-I. 



ds 



Zt = z, Xt=x 



e-KT-t) A\P{Xt))Z\T) 
1=1 

m n rn 

-re-'-(^-*) A\P{X,))Z\s) + e-''^^"*) ^ ^ Z\s)cuai{s, X,) 



i=l 1=1 



n „ m 



+ > I e 

k=i 



(3.16) 



v=l 



ds 



Zt = z, Xt= X 



By collecting the terms containing Z^ in sup„e and taking out the jth term in the supremum over 
u, we obtain 



dV{t,z,p) 
dzi 



-Ai{p) 

+ lim - jsupE \e-''^'^-*'> AUP{Xt))ZUT) 
e y I 



-r 

+ E / e-^^'-'^z-^{s)Y,iivP''{.i^''\s,x,-,e)>\de) 

/c=l"'* v=l 



ds 



Zt = z, Xt=x 



+ sup /i — sup /i 



-supE 

u 



e-^'^'^-i)A^{p{XT))Z'{T) 

n 

.^e-(«-*)^J (P(X,))Z^" (s) + e-'-(«-*) 5] Z^" (s)c,-,a, {s,X,) 



(3.17) 
where 
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I 1 -^IK ^. — 1 



E 



fc=l 



— z, Xt — X 



Y,A\P{Xt))Z\T) 



1=1 



1=1 



i=l 1=1 
1^3 



n „ rn 

-E / e-'-(-*)E^'(^) 
fe=i-^" 1=1 



1=1 



(3.18) 



Y,qi.p"'{i''\s,x,-,e)V{de) 



ds 



Zt = z, Xt = X 



Furthermore, define 



:= E 



e-'^(^-*)A^(p(XT))zj(r) 

7^ r n 

-re-'-(^-*)^^(P(X,))Z|(s) + e-'-(^-*)E^J(*)cj^"«(^'^^) 

ft „ 

+ E / e-'-(^-*)^|(.) 



i=l 



k=l ■ 



(3.19) 



and 



ds 



e-^(T-'^A^{P{XT)){Z^iT)+e) 

-re-'-(«-*)A-'"(P(X3))(Z^» + e) + g-'-t^-*) E(^'(^) + e)cj»a.(s, X.) 

+ E / ^"^'^"^ (^^■^)+^) 



k=l • 



(3.20) 



Y^qi.p^{i^''\s,x,-,e)y\de)]ds 



v=l 



Zt = z, Xt ^ X 



Then we can write (3.17) as 
dV{t,z,p) 



dzi 



A\p) 



14 
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1 

e-i-0 e 
i-T 



lim - \ sup/i - sup/i + sup/e - sup/o + eE e^'^C^"*) (P(Xt)) 

n 



i=l 



k=l 



-r{s-t) 



v = l 



(3.21) 



Since 'ZU we have that sup„e Ii — sup,^ /i < 0. By (3.151 there is an injective map between 



each functional 1^ and Iq for arbitrary such that ^ /q, hence sup„e /e — sup„ /q < 0. 
Consequently, 

dv{t,z,p) 



E 



v=l 



(3.22) 



X P''(7W(s,X,-,e''0)^'(de') 



Zt = z, Xt = X 



By applying the Ito formula to the process s i— ?> e '^'^ A' {P{Xg)) and taking conditional expectation 



with respect to Fs, we find that the right-hand side of (3.22) is zero 



□ 



In the case of a one-dimensional decision variable, i.e. an option of one commodity, this states 
an intuitively obvious result, namely that in the presence of an effective volume constraint, the 
usage of the option will lower its value. Note that if M > uT there is an e such that W — U and 



the map (3.15) is bijective. Hence, we obtain the result of Corollary 3.2 



4. The HJB-equation 



In the previous section, we studied the dynamic programming problem ( 2.5 ) first in the absence 



of an effective volume constraint for commodity i. We showed that in this case the optimal exercise 
rule can be determined explicitly and that the option does not lose value if used for this commodity. 
We also considered the problem in the presence of an effective volume constraint and showed that 
in this case it loses value when used. In this section, we determine an optimal exercise rule in 
the presence of an effective volume constraint. To this end, we first derive the associated HJB- 
equation. 
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For the reminder of the paper, we change the notation on the value function. Since rank{B) — 
TO, from now on we may write V exphcitly as a function of the factors X instead of the price 
P{x) = Bx, that is, we write V{t,z,x) instead of V{t,z,p) where the domain of V is modified 
accordingly. 



4.1. Necessary conditions. We derive now the HJB-equation of the problem (2.5). To this end, 
assume that value V exists. Then the Bellman principle of optimality yields 



(4.1) 



V(t, z, x) = sup E 

u<£U 



/•w ^ 

1=1 



Zt ^ z, Xt ^ X 



for all times < t < w < T. Rewrite the equation (4.1 ) as 



(4.2) 



sup E 



/ e-^' A'{P{X,))uids + e-™y(u;, Z^,X^) 
1=1 



-e-^'V{t,Zt,Xt 



Zt^z, Xt ^ X 



Furthermore, assume that V G C^'^'^(5). Then we obtain by the Ito formula 

/w 
d{e-^'V{s,Zs,X,) 

= / [e-'-^F,(s,Z„X,)-re-"l^(s,Z„X,) + e-'^^^y,,(s,Z„X,)a,(s,X,)]ds 

^=l 

e-"- Y,{aa^\jV^.^, {s, Z,, X,)d{W\ W^), 

+ / e-'''YV^i{s,Z,,X,)u[ds 
1=1 

+ / e-"^^K"(s,Z„X,)a,(s,X,)dM^, 



i=l 1 = 1 



k=l 



(4.3) +/ ^ / e-" [y(s,Z„X,+7W(s,X„e'=))-F(s,^.,X.)l A^'=(ds,d^'=) 



Here, aidWg = Yl,j ^ijdWl and 



where 7*^ is the jk'.ih element in the matrix 7. By compensating the Poissonian stochastic integral 



in (4.3), we find under suitable L^-assumptions on a and T4. , see [Ij, Thrm. 5.3.6, that the 
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Brownian and compensated Poissonian integrals in (4.31 are martingales. Thus the equation (4.2 1 
yields 



= sup E 

ueu 



1 " 

2 ^i^'^'^)^]^x-xii-'^^Zs,Xs)pij - rV{s,Zs,Xs) 
J2 / [v{s,Z,,X,+^^'^\s,X,,e))-V{s,Zs,X,)\,y'^{de) 

rn 

J2{A'{P{Xs)) + Vz^{s,Zs,X,))', 
1=1 

Define the integro-differential operator C on €^•^^'^{3) as 

n 1 ^ 

£F{t,z,x) = Ft{t,z,x) + ^ai{t,x)F^i{t,z,x) + - ^{aa'^)ijF^i^, {t, z, x)p.i^ 

i=l i,i 

(4.5) [F{s,z,x + j^''\s,x,e))-F{s,z,x)\i^''{de), 



(4.4) 



ds 



Zt = z, Xf = x 



k=l ' 



and rewrite (4.4) as 



= sup E 

ueu 



-r(s-t) 



\C~r)V{s,Z,,X,) 



m 

-J2{A'iP{Xs)) + V,,{s,Zs,Xs))' 



1=1 



ds 



Zt ^ z, Xt ^ X 



Under appropriate conditions on V, see, e.g. we can pass to the limit w t and obtain the 
HJB-equation 



(4.6) (£ - r) V{t, z, x) + sup |^(^'(P(x)) + V,. (i, z, x))u\t) | = 0, 

where the u varies over the set of IR™-valued functions defined on [0, T] satisfying the conditions 

0<u\t)<u\ [ u\s)^z\ [ u\s)ds<M\ 
Jo Jo 

for alU = 1, . . . , m and t e [0, T]. 



We observe from the equation (4.6 1 that the sign of quantity A''{P{x))+V2i{t,z,x), I — 1, . . . , m, 
determines whether the option should be exercised or not. From economic point of view, this 
quantity has a natural interpretation. Indeed, for a given commodity I, the function A''{P{-)) 
gives the instantaneous exercise payoff whereas the function V^i measures the marginal lost option 
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value. If the payoff dominates the lost option value for a given point (i, z, x) and commodity I, the 
option should exercised at the full rate. That is, for each commodity the option should exercised 
according to the rule 

if A'(P(a;)) < 

We also point out that this rule is in line with the case when there is no effective volume constraint. 
In this case, the marginal lost option value is zero and, consequently, the option is used every time 
it yields a positive payoff. In particular, we find that the presence of an effective volume constraint 
postpones the optimal exercise of the option for a given commodity I. 



4.2. Sufficient conditions. In this subsection we consider sufficient conditions for a given func- 



tion to coincide with the value function (2.5). These conditions are given by the following verifi- 
cation theorem. 



Theorem 4.1. Assume that a function F : S — >■ K satisfies the following conditions: 

(i) F{T,;-)^0, FeC'-'^HS), 

(ii) (£ - r)F{t,z,x) + Y.Zii'^^iPi^)) + Fz'{t, z, x))u[ < for all {t,z,x) € S and u e U, 



where C is defined in (4.5), 
(iii) The processes 

a) 9 ^ e— Er=i F.' Zs,X,)a,is, X,)dWs, 

b) ^ ELi L [^(*' Z,,X,+ 7^=) (5, e)) - F{s, Z,,Xs)] Nf^ids, d^), 
are martingales with respect to F. 

Then F dominates the value V . In addition, if there exist an admissible u such that 



{C — r)F(t, z, x) + sup 



Y,{A\pi^))+F.'{t,z,x))ui 



.1=1 



(4.7) 



(£ - r)F(i, z, x) + Y,{A\p{x)) + F,, (t, z, x))u[ = 0, 



1=1 



for all {t, z, x) G S, then u — u* and the function F coincides with the value V . 



Proof. Let u gU and t e [0, T]. By applying the Ito formula to the process t^ e ^^F{t, Zt, Xt), 



we find in the same way as in (4.3) that 



-''^V{T,ZT,XT)-e-'-'F{t,ZuXt) = d{e-^' F{s, Zs, X, 



)) 
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-'t i=i 

+ f X,)F,.,, (s, Z,,X,)d{W\ W^), 

„T rn 

z=i 



=1 1=1 



k=l 



By using the assumption (i), definition of the operator C and the equation (3.121, we obtain the 
equahty 



= 



^'F{t,ZuXt)+ I e-^'{L-r)F{s,Z,,Xs)ds+ j e"''^ ^ F,, (s, Z„ X.^ds 

Jt Jt 

nT n 

/ e-^'y2F,.{s,Zs,Xs)a,dW 
Jt 



k=l 



/ e-'-^ i^(5,Z„X,+7«(s,X„e'=))-F(s,Z„X,) N'^{ds,de). 



Conditioning up to time t and the assumption (iii) yields 



e-'-*F{t,Zt,Xt)+E 



^ e-'''{£-r)F{s,Zs,X,)ds 



Zt — z, Xt — X 



E 



1=1 



Zt = Xt = X 



By assumption (ii), we get 



(4.8) 0<e-'''F{t,Zt,Xt)--E 



Te— f;A'(P(X,))^ids 



Zi = z, Xt = X 



for all u £ U. Thus, the first claim follows. Now, if there exist an admissible u such that (4.7 1 



holds, then we would get equality in (|4.8|), i.e. for all w 

= e-'-*i^(i,Zt,Xt)-E 



/ e-"^A'(P(X,))^iids 
;=i 



Zt — z, Xt — X 



The conclusion follows. 



□ 
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5. Examples 

In this section we consider three examples. These examples illustrate two main issues. First, we 
compare two one-factor models (Example 1 and Example 2) with underlying Ornstein-Uhlenbeck 
factor dynamics, where the first model has a single Brownian driver whereas the other is driven 
by a sum of a Brownian motion and a compound Poisson process. To illustrate the effect of the 
jumps, the parameters of the factor dynamics are fixed such that the volatilities and the long 
term means are matched. In the third example, we study a two-factor model with underlying 
Ornstein-Uhlenbeck factor dynamics. As we will observe, the boundary conditions in the factor 
price dimensions are a delicate matter in this case. In all these examples, the aim is to find 



(5.1) y(t, z, x) = sup E 



Zt = z, Xt = X 



The boundary conditions in the x-direction are found by the same arguments as in [1]. The 
terminal condition is 



(5.2) 



V{T,z,x) = 0, 



for all z and x. This follows directly from the definition of the value function. Furthermore, the 
boundary condition in the z-direction, i.e. z — M, is 



(5.3) 



V{t, M, x) = 0, 



for all t and x. This follows from the fact that when z = M the only exercise rule available in U 



is the trivial one. The conditions (5.2) and (5.3) hold for all three examples below 



5.1. Example 1. Let the factor dynamics X be given by 

(5.4) dXs^K{fi~Xs)ds + adWs, Xt = x 

and P{x) = X, where s > t. Then it is well known that at time s > t, the solution 

(5.5) Xs ^ {x- ^i)e-''^'-'^ + ^i + a e-''^''-"^dWy, Xt^x. 
Furthermore, 

Xs r^N[^L + e-<^-'\x - m), ^(1 - e-2-(-*))^ . 
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With this specification, the value function (5.1 ) is given as a solution to the HJB-equation 



(5.6) Vt{t, z,x) + K{fi-x)Vx{t, z,x) + l-aVxx{t,z,x)-rV{t, z,x) + sup[{x-K + Vz{t, z,x))ut] = 0, 
with boundary conditions conditions in x-direction given by 

(5.7) V(t,z,x^ax)^uJ^ e-''^''-'^[{x„,a.- fi)e-''^'-'^ + fi]ds, 
and 

fT 

(5.8) V{t,z,x^,n)=u e-'-(^-*)[(a;™„-Ai)e-«(^-*)+M]ds. 

Je 

We remark that this is similar to the example in [J], Appendix A. However, we consider an 
arithmetic OU-process whereas in 0] the dynamics are given by an exponential OU-process. 



5.2. Example 2. To illustrate the effect of the jumps in the factor dynamics, we add in this 



example a compound Poisson process to the factor dynamics defined in (5.4) and match the 
expectation and volatility with Example 1. More precisely, consider the factor dynamic given by 
the Ito equation 

dXs = - Xs)ds + adWs + dY^, 

where the compound Poisson process Yg has Levy measure iy{dy) — fae^'^^l^yyojdy with /, a > 0. 
This equation can be written as 

(5.9) dXs = K{fi+ I yiy{dy)-Xs)ds + adWs+ [ yN{dy,ds), 

Jo Js. 

where the compensator 

/ 

yi^{dy) = -. 

a 



With this specification, the value function (5.1 ) is given as a solution to the HJB-equation 

1 f°° 

= Vt{t,z,x) + K{fL~ x)Vx{t,z,x) + -aVxx{t,z,x) + [V{t, z, x + y) ~ V{t, z, x)]fae "''"dy 

^ Jo 

(5.10) -rV{t, z, x) + sup[{x - K + V^{t, z, x))ut], 

ueu 

with boundary conditions in x-direction given by 

(5.11) V{t,z,xr,,a.)^u I e-'^(^-*)[(a;„,a,-/i)e-''(^-*)+/i+^(l-e-«(^-*))]ds, 

Jt Ka 
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and 

(5.12) V{t, z, xma.) =u r e-'-(^-*) [{x^a. ~ M)e-'=(^-*) + /i + ^(1 - e-'^(^-*))]ds. 



The solution to (5.9 1 is given by 

(5.13) X, = (x - - - At)e-''(''-') + M + - + a / e'^^'-'^UWy + f ( e-''^'-''^yN{dy,dv). 



It is easy to compute from the expression above that 



and 



nXs] = (x - ^ - /i)e-'=(^-*) + A + - , 
a a 



Var(X,) = E[(X, - E[X,])2] ^ ^(1 - e-^<^-'^){a^ + %] 



To match the volatility and the long term mean in Example 1 and Example 2, we solve the 
equations above for fl and a, when the expectation and variance is equal to that in Example 1. It 
follows that 

/ 

fi = fl , 

a 

and 



^ / 2 2/ 
a = \l (T^ ^. 

Then the mean and total volatility in Example 1 and Example 2 will be the same. This is good 
for comparison reasons, which will be discussed more in the next section. 

5.3. Example 3. The purpose of this example is illustrate the results when the price is driven 
by multiple factors. To this end, let Ai, /i^, /, k, K be non- negative constants. Consider the 

two-factor model 

dX{t) = a{t,X{t))dt + a{t,X{t))dW{t)+ [ j{t,X{t),ON{dt,d^), 

where dWt = {dWl,0), N{dt, d^) = (0, N^{dt, d^^^))) and the Levy measure is i^{dy) = (0, fKe-''yi{y>o}dy). 
Here, / is the jump frequency and k is the parameter of the exponentially distributed jumps. Fur- 
thermore, 

a{t,Xt) = {fi'-p'Xl-P^Xf), 

Ai 




Ht,x{t)) 
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e 



In component form we have, 



and 

dXl = - 

The solutions can be written as 



^^X^dv+ I ('N\dv,de), Xt=x'. 

R 



(5.14) 
and 

(5.15) Xl = e-'3'(^-*)x2 + j^^ j e-^"^'-^^eN^{de,dv) + 
To set up the valuation model, define the price function P : R'^ 



ev\dv,de 



as P{x) = a;^ + . Further- 



more, let Z be as in (2.3 1 with m = 1. The payoff is of call option type, i.e. A{p) = p — K. Then, 



the value function (2.5 1 reads as 



V{t, z, x) = sup E 
lieu 



-r{s-t) 



{P{Xt) - K)u,ds 



Zt = z, Xt= X 



From Proposition [3]l] we see that in the absence of an effective final volume constraint the optimal 
exercise policy u* is given by 

u iiP{Xt)>K, 
iiP{Xt)<K. 

for all t G [0, T], Hence, it is optimal to use the option whenever the swing yields a positive payoff. 
This is in line with 4J, in which no jumps are considered. 

Consider now the case with an effective volume constraint. The value function can be written 
in the component form as 



(5.16) V{t,z,x\x^) ^ supll 

lieu 



e-'-(^-*) [X] + Xl - K)usds 
This function is given as the solution to the HJB-equation 



Zt = z, Xf^x 



Vt{t,z,x\x^) + (^^ - l3^x^)V,^i{t,z,x\x^) - l3^x^V,^2{t,z,x\x^) 
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JR 

(5.17) + l-X'^V^ixi{t,z,x'^,x'^) - rV{t, z,x^,x^) + sup[(a;^ + x^ - K + V^{t, z,x'^ ,x^))ut\ = 0. 
2 usw 

With boundary conditions in t- and z-direction 



V{T,z,x\x^) and V{t,M,x\x^) =0, 



and boundary conditions in x-direction 



(5.18) 



(5.19) 



V{t, Z, Xjyl^^^, X ). 



To solve the problem (5.17l-(5.19l, we need to find the boundary conditions ( 5.18 1-( 5.19 1. We 



assume that we only have positive finite jumps, i.e. x^ > and that < << fi^. That is, 

the problem is solved in the first quadrant in the cc^x^-plane. 

Remark 5.1. The reason for choosing these spatial boundaries is due to the properties of the HJB- 
equation. In the x^ -direction we have diffusion, which requires boundary conditions at both ends. 
However, in the x'^ -direction we have transport in the positive direction, because the coefficient 
in front of V^^ is negative and that PDE is solved backward in time, and therefore no boundary 
condition is needed at x^ — x'^^^. At x^ = the derivative in x'^ -direction vanishes, thus no 
boundary condition is needed. 



In what follows, the calculations rely on the fact that the underlying factor dynamics are 



Ornstein-Uhlenbeck processes. By plugging in the processes and X'^ given by (5.14) and 



(5.15), respectively, into the value function (5.161 and rearranging the terms we obtain 



V{t,z,x^,x^) = sup<E 



+ ^(l-e-'^^(^-*)) + ^(l-e-^^(-))-K)..d. 



Zt — z, Xt — X 



(5.20) 



E 



+ r f e-^"^'-^k^N\de,dv))u,ds 

Jt JR ) 



Zt^ z, Xt^x 
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To compute V{t,z, x]^^^, a;^), we plug x^ = x\^^,^ into (5.20). Since the volatilities of the processes 
and are not state dependent, we can, by choosing x]^^^ sufficiently large, expect the 
trajectories of the process + X^ to be decreasing until the maturity since both the processes 
X^ and X"^ tend towards their long term means, /i^ and 0, respectively. Then it is optimal to 
start to exercise the option immediately with maximum rate until z = M since a^^j^^, + x'^ is much 
larger than the long time expectation. This argument holds since we only consider positive jumps, 
i.e. x]^^^ + x"^ > x]-,-^^^. Thus, if we start the process + X'^ in [xl^^^^x"^), we can define an 
optimal control as 

(5.21) Ms = wl,g[(^t+M^](s)- 

Then we get 

V(t-, z, a;^jj^, X ) — u e ^ ' ^2^,„^^e ^ ^ + x e ^ 



(5.22) 



+ ^2 (1 - e-^'(^-*)) + ^(1 - e-'^^(-*)) - k) ds. 
np'^ p^ J 



Here, we used the fact that the control u as defined in (5.21 ) is deterministic. This enables us to 
use the Fubini theorem and the martingale property for the second expectation in (5.20). 

On the contrary, when we start the process X^ at s^m; assumed to be sufficiently small, X^ 
will increase until maturity. It is thus tempting to wait as long as possible before we use the 
control, cf. equation (4.6) in [4J. However, since x'^ can be very large and the jump frequency is 
state independent, we are unable to draw this conclusion. To deal with this issue, we proceed as 
follows. First, we assume that the value function is continuous for all (.t^, x^) £ and both x]^^^ 
and x^^^ are finite. Using this we consider the deterministic part of the process X^ + X'^ starting 



at (x^j„,x^) at time t, that is, the first expectation in the value function (5.20). We observe that 
the integrand is a continuous function in time and will assume a maximum (and minimum) on 
[t,T]. Suppose it has its maximum at a time tmax- Furthermore, assume that 

The deterministic part will dominate the whole process at tj^ax- 

Then, due to continuity, the deterministic part will dominate the process on an interval (^1,^2) 
that contains tmax- We then choose a control defined as 



Us = 'ulse[ti,f2](s). 
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We substitute this into the expression (5.20). This is a deterministic control so, again, we can use 



the Fubini theorem and the martingale property to get rid of the conditional expectations. Define 



(5.23) J{ti,t2) -.^u I e 

'ti 



Then 



r{s-t) 



ds. 



(5.24) 



yit,z,xl^,^,x'^) = max J(ti,t2) 
ti,t2e[t,T] 



which is found by solving the (deterministic) maximization problem: 
rt2 



max u e 
ti,t2e[t,T] 



r{s-t) 



+ ^(l-e-^^(^-*))-if 



ds, 



subject to 



(5.25) 



u{t2 - h) < M - z. 



We solve the limits ti,t2 numerically by using elementary calculus methods. That is, to solve the 
maximization problem, define 



g{s,t) := e 



K 



This is the integrand in (5.231. By differentiating (5.23) with respect to ti and t2, we obtain the 



first order necessary conditions 

(5.26) 

and 

(5.27) 



dt2 



dh 



ugit2,t) = 0, 



-ug{ti,t) = 0. 



Furthermore, at the boundary where t2 ~ ti + '^^^ ^ , we find 

dJ(ti,t2) , , M - z ^ 
(5.28) y^'^ " = u{g{h + -^,t)~g{h,t)) = 0. 
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To conclude, we solve ti and ^2 from these three equations, substitute these into (5.231 and see 
which gives the highest value still satisfying the constraint (5.25). 



6. Numerical experiments 

In this section we will present the numerical solutions of the HJB-equations from the three 
examples in Section 5. All equations have the boundary conditions that the option value is zero 
when t = T and z — M. In addition we truncate the boundary in infinity, and the truncated 



boundaries requires boundary conditions. In Example 1 we solve (5.6) with the boundary con 



dition (5.7) and (5.8) on the truncated boundary. Similarly in Example 2 we solve (5.10) with 



the boundary condition (5.11) and (5.12) on the truncated boundary. And finally in Example 3 



we solve (5.171 with the boundary condition (5.22) and (5.24 1 on the truncated boundary. Be 



low we specify the model parameters, which are chosen for the purpose of illustration, and the 
discretization parameters of the numerical scheme. 

6.1. Numerical scheme. In Example 1 and Example 2 the HJB equations are PDEs defined 
over the variables t, z and a;, while the HJB equation in Example 3 is defined over the variables 

z, and . The equations are solved with finite difference methods (EDM). We use a first 
order Euler scheme in t-direction. The z- and x^-directions are handled explicitly with first order 
upwind schemes, while the x or x^-directions are handled implicitly with a second order central 
difference scheme. We start the time stepping at t = T and go backward in time until t = 0. 

The domain is discretized with a uniform grid in the i-, x^- and z-directions whereas in the 
X- or a;^-direction we use an adaptive grid. The integral term is approximated with numerical 
integration, more precisely, the rectangle method with second order midpoint approximations. 
The truncated boundary in the direction of jump, i.e. the x-direction in Example 2 and the x^- 
direction in Example 3, causes some problems for the approximation of the integral, which is 
supposed to have upper limits at infinity. This problem is solved by linearly extrapolating the 
option price outside the truncated domain, and integrating up to a level where we get sufficiently 
accurate approximation of the integral. 



6.2. Numerical examples. In this subsection we study the numerical examples from the previous 
section. They are all motivated by some swing options traded in the Scandinavian electricity 
market, which are called " Brukstidskontrakt" . Such a contract gives the owner the right to buy 
a certain amount of electricity for her own selection of hours during 1 year. More precisely, for 
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each of the 8760 hours in 1 year the holder of the contract must choose whether or not to use the 
contract. In our example we set the portion to 50%, i.e. the holder must choose 4380 hours. 

In practice the contracts are usually paid in advance and not for each time it is exercised, so 
the strike price will he K — in all the examples. We also use T — 1, M — ^ and u — 1. The 
three examples are further specified in the following. 

Example 1: The model is specified by k = 0.014, /i = 40, cr ~ 2.36. The truncated domain 
is defined by Xmin — ^8.7,Xmax — 61.3. The discretization parameters are At — Az — The 
grid in x— direction is adaptive and consists of 671 grid points, with higher grid point density 
around fi and lower density near the truncated boundaries, Xmin and Xmax- 

Example 2: The model is specified by k = 0.014, a = 0.4, / = 0.04, ^ = 39.9, cr = 2.3387. 
With these parameters the model i Example 1 and Example 2 have thesame mean and volatility. 
The grid and the truncated domain is as in Example 2. 

Example 3: The model is specified by ^ 0.014, ^ = 40, cr = 2.36, /J^ ^ 0.04, k ^ 0.014, / = 
0.04. The truncated domain is defined by a;^„j„ = 17.2, a;^^^ = 62.8, = 0,x^^^ = 9. The 

2 

discretization parameters are At — Az = gj^, Ax^ — . In the a;^-direction we use an 
adaptive of 1200 grid points, with higher grid point density around fi and lower density near the 
truncated boundary in x^-direction, i.e. — and x^ — xl^^,j.. On the truncated boundary 



condition in x^-direction we need to calculate (5.22) and solve (5.24). With the parameters in our 



example, ti and t2 in (5.24) turn out to he ti — t and t2 = t+ — z. The truncated boundary in 
x^-direction needs no boundary condition due to the nature of the PDE. We also tried adaptive 
grid in the x^-direction, but this did not seem to improve the accuracy. 

In the following we visualize the numerical solution of these three examples. The two things 
we are most interested in are the option prices and the trigger prices. The trigger prices are also 
referred to as exercise curves, and they tell us when to exercise and when to hold. The option 
price is a function of two variables for each point in time in Example 1 and Example 2, and can 
then be visualized in a 3D plot. However, the option price in Example 3 is a function of three 
variables and is therefore harder to visualize, even for a fixed point in time, but we present it in 3 
plots. The trigger prices are in [3j presented as exercise curves (see more in figures below). These 
prices are presented similarly for the results of Example 1 and Example 2. But for Example 3 
the trigger price is actually a 3D surface for each point in time. We solve this by projecting the 
surface down to two different planes. We could also have made 3D surface plots of the trigger 
price in this example, but we think 2D plots are more instructive. 
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Figure [T(a) and Figure 1(b) show the option price at t — 0.5 for Example 1 and Example 2 
respectively. We see that the two plots a quite similar, and that the option price increases with 
increasing x-values. This makes sense from an economical point of view, since one would expect 
that a higher spot price results in a higher option price. Mathematically we see it from the value 



function (5.1 1 and the fact that the solution functions (5.51 and (5.131 are increasing functions of 



X. Furthermore, the option price decreases with increasing z-values which is in accordance with 
Proposition |3.3[ stating that whenever using the option it loses value, which makes economical 
sense. 



Figure 1(c) shows the difference in option price between Example 1 and Example 2 at the 



same time level. More precisely Figure 1(c) shows the price from Example 2 minus the price 



from Example 1. We see that this difference is negative, which means that the option price is a 
little higher when we assume an underlying jump process. This is reasonable since the Gaussian 
OU-process in Example 1 has a symmetric distribution whereas the non-Gaussian OU-process in 
Example 2 has a positively skewed distribution. This positive skewness, which increases value in 
financial markets, is caused by the fact that we only have positive jumps. 





(a) Option price Example 1 



(b) Option price Example 2 



T^l Jit 1 H L 1 E J 111 I jtiJ E H [ Ih 




(c) Difference between option prices in Example 
1 and Example 2 



Figure 1. Option prices in Example 1 and Example 2 
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In Figure [2] we show the exercise curves for Example 1 and Example 2 at time = 0.5. The red 
curve corresponds to Example 2 and the black curve corresponds to Example 1. We see that the 
red curve lies more to the right than the black curve. 
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Figure 2. Exercise curves for Example 1 and Example 2 at time = 0.5. Exercise 
when at points to the right of the curve and hold when at points to the left of the 



In Figure 3(a) - 3(c)| we show the option price of Example 3 at i = 0.5. The function is plotted 
as a function of a:^ and z for three values of x'^. We see that for each value of the plot looks 
similar to the plots in Figure 1(a) and Figure l(b)[ only the level of the surfaces changes a little. 



In Figure 4(a) - 4(b) we visualize the trigger price for t = 0.5 in Example 3. Figure 4(a) shows 
the exercise surface projected down to the price, z-plane (where price — +a;^), for various values 



of X . Figure 4(b) shows the exercise surface projected down to the a; x -plane, for various values 



of z. 



It is worth noting that the slope of the curves in Figure 4(b) is approximately — = —0.35. 



For example if we study the red curve and remove the point where x^ — and make a linear least 
squares approximation of it, it will have a slope of —0.3517. This means that the negative ratio of 
the mean reversion speeds approximates the slope of the exercise curves in the x^x^-plane. This 
is plausible from an economical point of view for the following reason. If the mean reversion speed 
is smaller for x^ than for x^, the holder will exploit a deviation from the long term mean earlier 
for process X"^ than by exercising the option. That is, she would require a higher contribution 
from x^ to the price than from x^ before exercising. This is because a high value of x^ is likely to 
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reduce more quickly and therefore it is beneficial to exercise with a lower value of in relation 
to . On the contrary, if has a high price it is more likely to stay high longer. In this case, the 
holder might wait for even higher prices. A similar reasoning can be done for the opposite case 
when the mean reversion speed is bigger for x^ than for x^ . 

Notice that at t = 0.5 and for any value of z, the holder should always hold if the values of 
x^ and are small enough. The reason for this is the following: if the underlying price is small 
enough, you would expect it to be higher than this the rest of the time, and for t = 0.5 it will 
therefore be beneficial to hold since M — 0.5. However, in Figure 4(a)| it can be seen that the 
exercise curves seems to stay above the line z=0.G01. This is due to numerical error/instability. 
For finer grids the exercise curves will be closer to z = at for small values of P ~ + X'^ and 
t — 0.5. Similar effects can be observed for other values of t. The solution to this problem is either 
higher resolution on the grid, which requires more memory on the computer, or more accurate 
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numerical schemes. The curves in Figure [4(b) [ bends a Uttle in the lower right corner, and this is 
due to the mentioned instability in the numerical scheme. 
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(a) Exercise curves plotted in the price, z-plane (b) Exercise curves plotted in the x^x^-plane 

Figure 4. Trigger prices in Example 3 



6.3. Numerical accuracy in Example 3. The analysis of the numerical scheme in Example 1 
and Example 2 is similar to that of [4 . To analyse Example 3 we study the Courant-Friedrichs- 
Lewy (CFL) condition. The CFL number of the HJB equation in Example 3 is 

1 1 

_ 3200 n 0/1^2 , 3200 1 



40 3198 

= 320o(l4+^l^«, 

A necessary condition for convergence is that the CFL number C < Cmax, and in our case with 
e (0,9) we see that the CFL number C < 1. As mentioned we have observed some small 
instabilities in the numerical solution. We have also tried with larger values of Az compared to 
At, but the reported discretization parameters seem to give most accurate solutions. In order to 
establish convergence an implicit scheme should be developed, but this is not done in this work. 

In the following we will present some evidence that the numerical solution converges to the 
correct solution of the HJB equation. We attempt to evaluate both the numerical scheme and the 



calculated boundary conditions using (5.19). This is done by trying to see how well they fit for 
extreme values of x^, i.e. x^ >> fi. We have used a very fine grid to solve HJB equation where V 
is required to be linear in x^-direction at the boundaries. This is similar to :4!, where it is shown 
that this type of inaccurate boundary condition gives quite accurate solutions. Now this solution 



can be compared to the values we get from Equation (5.22 1. 
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The difference between the numerical solution of the HJB equation and the value calculated 
by (5.22) is illustrated in Figure [S] We see that the difference is relatively small compared to the 
option value. For this example with t = 0.5 and z = 0.4 it is about 6 orders of magnitude lower 
than the option value. This is sufficiently small for us to trust the numerical solver. The difference 
may come from all of the following five sources: assumptions that the control is as described in 



(5.21), numerical inaccuracy/instability of the scheme, truncation in direction, truncation in 
direction, extrapolation in x'^ direction and linear boundary condition on the PDF solver. 



Difference between numerical solution and boundary condition at time = 0.5 and i = 0.4. 




Figure 5. Difference between the numerical solution and (5.19) for large values of x^ 



7. Conclusions 

In this paper, we developed and analyzed a valuation model for swing options on multi- 
commodity markets. The model is formulated as a dynamic programming problem, where the 
underlying dynamic structure is given by a multi-dimensional Levy diffusion. This process models 
the price evolution of the commodities. The commodity prices are driven by a multi-dimensional 
Brownian motion and a multi-dimensional compound Poisson process. The introduction of the 
compound Poisson process is important since it allows non-Gaussian price evolution. This is im- 
portant, in particular, on electricity markets, see, e.g. [31[T31[T7]. Furthermore, this model allows 
us to take into account jumps in price processes, which is also important on electricity markets. 

From a analytical point of view, this study provides a multi-dimensional generalization of the 
analysis in |4j. First, we analyze the model in the absence of an effective volume constraint. Along 



SWING OPTIONS IN COMMODITY MARKETS 33 

the lines of [1], we find that in this case the option does not loose value if used. Moreover, we 
prove that in the presence of an effective volume constraint for a given commodity, the usage of 
the option for this commodity will lower the value of the option. This is a intuitively appealing 
from the economical point of view. To tackle the problem of finding an optimal exercise rule and 
the price of the option, we analyze the pricing problem using the Bellman principle of optimality 
and derive the associated HJB-equation. In Section 4, we obtained an optimal exercise rule which 
states that if the immediate exercise payoff dominated the lost option value for a given commodity, 
then the option on this commodity should be exercised at a full rate. In particular, we conclude 
that this optimal exercise rule is a bang-bang rule. We also provide a verification theorem, which 
states conditions under which a given function coincides with the value function. 

In addition we illustrate the results with three examples which we study numerically. We set up 
a straightforward FDM scheme to solve the associated HJB-equations. The numerical experiments 
seems to give reasonable results from both mathematical and economical points of view. In the 
last of our examples we have also given evidence for convergence of the numerical solution. 
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